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Abstract 

An effective interaction between fermions in a Bose-Fermi mixture is derived. It is induced by 
density fluctuations of the bosonic background. The contributions from states containing both 
one and two virtual phonons are taken into account self-consistently. The time dependence of 
the effective interaction has been removed by assuming that the velocity of the fermions at the 
Fermi surface is much larger than the sound velocity in the Bose gas. This assumption is more 
appropriate for the actual experimental situations than the usual approximation of neglecting 
retardation effects. The effective interaction turns out to be attractive and, as a consequence, can 
give rise to a superconducting phase in the Fermi component of the mixture. The fermions are 
considered in only one magnetic state, so that pairing can be effective only in odd-Z channels. It 
has been found that the onset of the superconducting phase can occur at temperatures ( > 100 nK ) 
of the same order of magnitude as the Fermi temperature ( ~ 300 nK ), and the energy gap in the 
excitation spectrum is a small fraction ( ~ 1% ) of the Fermi energy. 
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I. INTRODUCTION 



Effective fermion-fermion interactions induced by the exchange of virtual phonons in 
fermion-boson mixtures were already studied in the 1960s, in connection with dilute solu- 
tions of 3 He in superfluid 4 He, see e.g. Ref. Recently, a renewed interest in this subject 
has arisen, following the actual availability of trapped atomic gases at a temperature suffi- 
ciently low to exhibit quantum properties. In particular, physical situations are reproduced, 
in which a highly condensed Bose gas is mixed with an almost degenerate Fermi gas and 
the interaction between these two gases can be studied at various values of temperature and 
density 0, 0, HI 11 • On the theoretical side, particular attention has been devoted to the 
possibility of the onset of a supercoducting phase in the Fermi component of the mixture 
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For the fermion gases used in experiments the bare interatomic potential is 
repulsive. Then, the necessary attractive interaction between two fermions should be pro- 
vided by the exchange of virtual phonons, through a mechanism similar to that responsible 
for the superconductivity phenomenon in metals. Investigations of the physical conditions 
for a a Cooper instability in the s-wave and p-wave channels have recently been performed 

aaai. 

For a bare fermion-boson interaction with zero-range, the induced fermion-fermion in- 
teraction is essentially proportional to the correlator of density fluctuations of the bosonic 
background. In the references [H H, H> El> the density fluctuations have been calcu- 
lated taking into account only single- phonon intermediate states. In the present paper the 
contribution from two-phonon intermediate states is also included. This does not simply 
amount to an improvement of the approximation scheme, because two-phonon and one- 
phonon states occur on the same basis in the evaluation of the effective fermion-fermion 
interaction, at lowest order in the fermion-boson coupling constant. Moreover, when the 
fractions of condensed and non-condensed components of the boson gas are comparable, the 
one-phonon and two-phonon contributions turn out to be of the same order of magnitude. 

The induced interaction between fermions is given by the boson density-density correla- 
tor, evaluated for the time interval between two successive fermion-boson interactions. In 
order to simplify calculations, different approximations can be introduced to eliminate this 
time dependence. In Refs. (H 0, 0, 0, retardation effects have been neglected, always by 
setting the frequency to zero in the time-Fourier transform of the correlator. For a degener- 
ate Fermi gas, this approximation can be justified if the Fermi velocity is much lower than 
the sound velocity in the Bose gas. In the present approach the opposite approximation 
is used: the correlator is taken at equal times, i.e. the induced interaction between two 
fermions is mediated by the static correlation function for boson density fluctuations. This 
approximation is more appropriate for the physical conditions obtained in the mixtures of 
trapped Bose-Fermi atoms produced in the laboratory @ 5 S S ■ F° r these mixtures the 
density of the degenerate Fermi gas is a relevant fraction of the density of the Bose gas. 

The aim of the present work is to assess the possible occurrence of a superconducting 
phase in the Fermi component of an ultracold mixture of Bose-Fermi gases, for values of 
density and temperature obtained in actual experiments, and to give an estimate of the gap 
in the spectrum of excitation energies. Although all the experiments have been performed 
with harmonic-oscillator traps, a spatially homogeneous system is considered. In this case, 
calculations are considerably simplified and compact analytical expressions can be derived. 

At the densities considered in the experiments, for the two species of gases we can assume 
that the interatomic interactions have a short range and mainly occur in the s-wave channel. 
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Moreover, according to the experimental situations, both the Fermi and Bose atoms are 
considered in only one internal state. In this case, the Pauli principle allows us to neglect the 
bare interaction between fermions, at least in a mean-field treatment. In fact, in the zero- 
range limit for the interaction, the contributions from the direct ( Hartree ) and exchange 
( Fock ) terms of the mean-field cancel each-other exactly, while short range correlations 
should not play an important role for a sufficiently dilute Fermi gas, i.e. for a FF np << 1 
where a pp is the scattering length and tip is the fermion density. Concerning the Cooper 
condensation, it can occur only for odd values of the orbital angular momentum for a spin- 
aligned couple of fermions, so that there is no contribution from a s-wave interaction. 

Experiments have been performed with gases, which show repulsive boson-boson inter- 
actions |H IH , this ensures the stability of the mixture towards separation of the two species 
of gases allowing us to study mixtures with a sufficiently high number of atoms. Moreover, 
the interaction between a fermion and a boson in the mixture 40 K- 87 Rb of the Florence 
experiment [H| is strongly attractive, thus the induced fermion-fermion interaction can be 
very effective to bind fermions in correlated pairs. The calculations are performed on a 
general ground, however, when more phenomenological aspects are concerned, they mostly 
refer to the physical environment of Ref. [^] . 



II. FORMALISM 

In this section the effective interaction between two fermions, mediated by density fluctu- 
ations of the Bose gas, and the equation for the energy gap in the excitation spectrum of the 
Fermi gas are derived. The calculations are performed by means of functional methods thor- 
oughly, since both the properties of the Bose condensate and the induced fermion-fermion 
interaction can be directly obtained in the same scheme ( for a introduction to functional 



methods see e.g. Ref. [10|, [ll| 



A. Effective interaction 

The starting point is the functional-integral expression for the grand-canonical partition 
function of the Bose-Fermi mixture, that reads 



V{<f>, 0*, V, r) exp [ - (S B {<f>, 0*) + S F (tf>, r) + Si(<P, 0*, V, P 



The integration is performed over the complex field 0(x, r) ( bosons ) and the Grassman field 
ijj(x, t) (fermions), which are respectively periodic and anti-periodic in the imaginary-time 
interval (0, (3 = 1/T) ( units such that h = c = k B — 1 are used ). Since we deal with spin- 
polarized particles, 0(x, r) and "0( x ; T ) are one-component fields. The action functional has 
been written as the sum of three terms: the action for the boson field 



S 



B{ 



dr \ dx 



^(x,r)( 
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(2) 



the action for the fermion field 



Sp(ip,ip*) = J dr i <7x 



r (x,r — 







dr Imp 



(3) 
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and the interaction term containing both fields 



= A f dr I rfx ^(x,r)^(x,r)0*(x,r)0(x,r) . ( I) 



o 



Here /ig and fip denote the boson and fermion chemical potentials respectively. In writing 
the action functional, the usual zero-range approximation for the interactions between atoms 
has been adopted The corresponding coupling constants are defined as follows: 

7 = iiiaBB /tub and A = lixasF l^fiR-, where cibb and cibf are the boson-boson and fermion- 
boson scattering lengths respectively, and trr = uibtuf / (ms + is the reduced mass for 
a boson of mass tub and a fermion of mass mj?. Finally, according to the previous remarks, 
in Eq. (jSJ) the fermion-fermion interaction has been omitted. 

The Bose gas is treated within the Gaussian approximation, i.e. all terms of order higher 
than second in the fluctuations of the boson field about the uniform and constant condensate 
amplitude are neglected. This corresponds to the Bogoliubov approximation for a simple 
Bose gas. For a weakly interacting gas, n B a BB « 1 ( n B is the boson density ) and at low 
temperatures, this approximation is well justified |l2|. The first step consists in evaluating 
the condensate amplitude about which the perturbative expansion can be performed. The 
standard methods employed for a simple Bose gas 0, 11 1 are extended to a fermion-boson 



mixture. First the action functional is integrated over the fermion fields, then, a uniform 
and constant boson field 0(x, r) = ip, which makes the action stationary, is searched for. 
The integration over the fermion fields can be performed explicitely, since these fields occur 
with a Gaussian form in the functional integral. This gives the two following contributions 
to the action functional 

-Tr lnGg^x-xV-r') - Tr In (l + AG x (x - x', r - r')0(x, r)0*(x, r)) , (5) 

where the trace operation implies integration over the continuous variables x, r. Here 
G'o 1 (x — x',r — t') is the inverse of the Green function for free fermions. The first vari- 
ation of the action functional takes the contribution 

—A 0(x, r) G(x, x, r, r + + ) = A 0(x, t)uf, 

from the second term of Eq. (pj. The quantity G(x, x', r, r') is the Green function for 
fermions interacting with the scalar field 0(x, r). Adding the contributions coming from the 
boson term of the action, for the density of the zero-momentum Bose condensate, n° B = p 2 , 
the following equation is obtained 

o — An F 
n B = • (6) 
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This equation has been derived assuming, without loss of generality, that the quantity ip is 
real. We can see that the effect of the interaction with the Fermi gas simply amounts to 
replacing the chemical potential /xg with the effective value \ib ~ Anp. 

The next step is the introduction of the change of variables 0' = <p — (p in the terms 
Sb{,<P,4>*) and Si(<p,<p*,ip,ip*) of the action functional. The term with p 2 represents the 
interaction of the fermions with the condensate mean-field, \n% and can be added to the 
chemical potential \xp in Sp(ip,ip*). The terms containing products of more than two am- 
plitudes <f)', are neglected according to the Bogoliubov approximation. Thus, the action for 
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the fluctuations of the boson field is given by 



1 rfi r * 

- / dr dr' / dxrfx' 0^(x, r).D _1 (x, x', r, r')<//(x', r') 

2 Jo 7 

+ ^ jf rfr y dx^'^x, r)f(x, r) + c.c. 



(7) 



where the vector field 



the source term 



y*(x,r) 



T(x,r) 



V>*(x, r)^(x,r) 



^(x,r)^(x,r) ; ' 
and the inverse of the Green function for the fluctuations <ft' 

D -1 (x, x', r, t') = D^fa x , r, r') + t/(x, x, r, r 

are introduced with the inverse of the Bogoliubov propagator 



/ d 



D n l (x, x',r, r'l 



V 2 



<9r 2ms 



\ 



\ 







dr 2m B 



+ jn° B 



5(t - r')5(x - x') 



and the coupling matrix 

^(x,x',r,r') = A 



^(x,t)V(x,t) 

^*(x,t)VKx,t) 



<5(t — r')5(x — x') . 



In order to obtain an effective interaction for the fermions, the integration over the fluctu- 
ations 0' and 0'* should be performed. The result is an effective action containing only the 
fermionic fields if) and ip* 



S eff (iJ),iJ)*) 



dr I dx 



\ dr 2mp 



{Hf - Ar^)Wx,r 



+1 b^detl^ 1 ) + l - In [det(l + £>„?)] 



1 fP 
2 



jf drdr'y" dxrfx'f t (x,r)Z?(x,x , ,r,r , )f(x / ,r') 



(8) 



The term b^det/J^ 1 ) together with the action evaluated at the stationary point <fi' = 0, 
give the grand potential in the Bogoliubov approximation for a Bose gas with the effective 
chemical potential \i% — \n F . 

The last two terms on right hand side of Eq. (JHJ) represent an effective interaction among 
fermions. When expanded in powers of the fermion fields, both terms give rise to n-body 
forces. Only terms containing up to four fields are taken into account, so both states 
with more than two virtual phonons in the boson-density fluctuations and n-body forces 
among fermions, with n greater than two, are neglected. We expect this to be a reasonable 
approximation for sufficiently dilute gases. According to this approximation, the boson 
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Green function is replaced with the Bogoliubov propagator, D = Dq in the source term, and 
terms up to second order are retained in the expansion 




(9) 



The space-Fourier transform of the contribution from the source term of Eq. (JHJ) to the 
effective interaction between fermions is given by 



Vifa, n - r 2 ) = -X 2 n°M n \q, n - r 2 ) + D^ n \q, r 2 - r x ) + 2D 12) (g, n - r 2 )\ , (10) 



where Dq are the elements of the 2x2 matrix Do; for their explicit expressions see, e.g., 
Ref. 0| ( the definition of D used here differs from that of Ref. 0] by an overall minus 
sign ). Concerning the contributions from the expansion of Eq. Q, the term with n — 1 
represents the mean field acting on the fermions due to the noncondensate fraction of the 
Bose gas. This term can be added to the chemical potential of the fermions, giving together 
with the analogous term of Eq. (JHJ) the effective chemical potential p,* F = ftp — Xn B for the 
fermions. The space-Fourier transform of the contribution from the term with n = 2 reads 



This term is the contribution to the effective interaction from states containing two virtual 
phonons. 

A rough estimate of the terms of Eq. @ corresponding to the exchange of n > 2 
virtual phonons between fermions can be done by means of the following considerations. 
The expansion of Eq. (JHJ) contains correlations of the boson density fluctuations of higher 
and higher order when n increases. According to the basic assumption made in this paper, 
the correlations will be evaluated at equal times. In this case, we should expect that the 
contribution of a term of the series with respect to that of the previous term, contains an 
extra power of Xn' B . Here, n' B is the density of the noncondensate bosons. In addition, 
when n increases by a unit, the effective action of Eq. (JHJ) contains an additional integration 
on the time, which gives a multiplicative factor of the order of magnitude as (3. Thus, we 
can interpret the right hand side of Eq. Q as a perturbative expansion in powers of the 
adimensional parameter \n' B /T. Similar considerations can be done for the terms neglected 
in the perturbative expansion of the last term of Eq. (JUJ). At the densities and temperatures 
considered in this paper the parameter \Xn' B /T\ takes a value of ~ 0.15. Since the order 
parameter of the superconducting phase, in the approach of this paper turns out to be 
not critically dependent on the fermion-fermion interaction ( see Eq. (J2"H|) of Sec. Ill ), 
corrections to the effective fermion-fermion interaction coming from terms of order higher 
than second will be neglected in order to simplify calculations. 

B. Pairing field 

The effective interaction obtained involves fermions at two different times. The approxi- 
mations which can be introduced in order to get rid of the retardation effects and simplify 




dk 



W i; (*,ri-r 2 W i; (k-q,r a -ri) 
+D { 12 \k,T 1 -r 2 )D { 12 \k- q ,r 2 -r 1 )]. 



(2tt) 3 



(11) 
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calculations, will be discussed later. Now, the equations for the effective field representing 
the pairing of two fermions, are derived. 

With the approximations discussed above, the partition function of Eq. (0) assumes the 
following factorized form 



Z = exp(-PQ B ) J V(i/;,r)eM-Seff(i[>,r)} 



(12) 



where Qb is the grand potential of the Bose gas, while the fermion effective action of Eq. 
dHl becomes 



Seff(1>,r)= I dT f dk[PQt,T)(— + — -/4)V(k,7 
1 rP r 

— dT X dT 2 J Y[ dki S(k 1 + k 2 - k 3 - k 4 



2(2vr) 3 Jo 

x K//(k 3 - k ls n - r 2 )V*(ki, ri)V*(k 2 , r 2 )</>(k 4 , r 2 )^(k 3 , r x )] , (13) 

with V e// (k 3 - ki,n - r 2 ) = Vi(|k 3 - k^n - r 2 ) + ^ 2 (|k 3 - k a |,n - r 2 ). 

The pairing field is identified with the composite field which makes the effective action 
S e ff(ip,ip*) stationary. This can be determined by means of a Stratonovich-Hubbard trans- 
formation 0, 14]| introducing auxiliary fields for the couples ipip and ip*ip*, and by equating 
to zero the first variation of the resulting effective action with respect to the auxiliary fields. 
After the transformation, the exponential in the integral of Eq. ()12j) is replaced by 



exp[ - S eff (ip,ip*)} 



cxp 



x I £>(A*,A)exp[^n dr i J II dk i 



1 r 



A*(ki, k 2 , n, r 2 )V eff (k 1 , k 2 , n, r 2 ; k 3 , k 4 , r 3 , r 4 )A(k 3 , k 4 , r 3 , r 4 ; 



+ -y (iki(ik 2 y dridr 2 [A*(ki,k 2 ,ri,r 2 )^(ki,r 1 )^(k 2 ,r 2 ) 
+^*(ki, r!)V*(k 2 , r 2 )A(k 1 , k 2 , n, r 2 ) 



(14) 



where A and A* are bosonic auxiliary fields and the matrix 

K//(ki, k 2 , n, r 2 ; k 3 , k 4 , r 3 , r 4 ) 
is the inverse of the effective interaction 

< kik 2 |V e// (ri,r 2 ;r 3 ,r 4 )|k 3 k 4 >= 



1 



(2tt) 



:K//(k 3 - ki, n - r 2 )<5(r 3 - ti)<5(t 4 - r 2 )5(ki + k 2 - k 3 - k 4 ) . 



In the functionals containing the coupling between fermions and auxiliary fields, only the 
parts of A and A* which are antisymmetric under the permutation 1^2 survive. This 
reflects the prescription of the Pauli principle for fermions in the same spin state. 
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The integration over the fermionic fields can be formally performed, giving the following 
action for the auxiliary fields 



S(A,A*) = - / J] dr { /n dk il[ A*(k 1 ,k 2 ,r 1 ,r 2 )K7K k i' k 2' r i' r 2; k 3 ) k4,r3,r 4 ) 

xA(k 3 ,k 4 , r 3 ,r 4 )] -iTrln(l-^ A), (15) 

where Qq is a 2 x 2 diagonal matrix, whose elements are the propagator for independent 
fermions with the effective chemical potential fi* F , and the matrix A is given by 



A(ki,k 2 ,Ti,T 2 ) 



A(k 1; k 2 , n,r 2 ) 

A*(ki,k 2 ,ri,r 2 ) 



The equation for the pairing field can be obtained by equating to zero the functional 
derivative of S(A, A*) with respect to either A* or A, equivalently. It reads 



A(k, k', t-t') = J dk x dk 2 K//(ki - k, r - r')6(k + k' - k x - k 2 ) 

x^ 12 )(k 2 ,k 1 ,r-r'), (16) 



where £ (12) (k 2 , ki,r — t') is the anomalous propagator for fermions interacting with the 
pairing field A ( see, e.g., Ref. ). Because of time-translation invariance both the 

pairing field and the fermion propagator depend only on the difference of times. 

It is convenient to introduce the center of mass momentum P = k+k' = ki + k 2 , thus, for 
a given value of P, an equation containing only one integration variable, k, can be obtained. 
The solutions of Eq. (|16|) with P ^ 0, correspond to a breaking of the space-translational 
symmetry, i.e. the LOFF phase [15j. Here, only solutions with vanishing center of mass 
momentum are considered. In this case, the pairing field and the anomalous propagator 
depend on only one momentum, k = — k' and ki = — k 2 , respectively. The propagator 
Q^{ki, t — t') is given by the Fourier series 

G^Hki r r>)= W e-^(^)A(k 1 ,o; w ) 

y 1 ' } ~ (3 2 ^^ n + e kl +A(k 1 ,u n )A*(k 1 ,co n )' [t) 

where u n are the Matsubara frequencies for fermions and £j? = (kl/2m F — [i* F ) 2 . 

Now, two different paths of approximation which allow us to neglect the time (or fre- 
quency) dependence of the pairing field, are discussed. Generally, for evaluating the time 
behavior of a propagator we should make an analytical continuation to real times ( —it — > t ) 
or, equivalently to real frequencies ( iuo n — > u ). However, for the two cases that will be en- 
visaged, we can directly refer to the equations just derived. Both approximations are based 
on the comparison of the characteristic times of propagation for fermions and for bosons. 
For a quasi degenerate Fermi gas only the particles with momenta lying about the Fermi 
surface partecipate in the pairing process. Then, for a given momentum q, the fermion 
time-scale can be estimated as t F = l/v F q. While for bosons, the quantity t B = l/c s q, 
where c s is the sound velocity in the Bose gas, can be assumed to give an estimate of the 
characteristic time. The comparison between the two time-scales reduces to compare the 
velocities Vp and c s . The case, when c s >> Vp, is examined first. This is the assumption 



usually made when studying boson-induced interactions between fermions for Bose-Fermi 
mixtures. For c s >> Vp retardation effects on the induced interaction can be neglected, i.e. 
only the u = component of the time-Fourier transform of the boson propagators in Eqs. 
()1U|) and (fTTj) is retained. Also, the fermion propagator (^ 12 )(k, r — r') can be evaluated at 
equal times (r = r'). The pairing field in turn becomes an instantaneous function of the 
time difference, A(k, r — r') = A(k)5(r — t'). The equation obeyed by the pairing field is 
the usual gap equation of the BCS theory. The second approximation is the opposite of 
that just discussed: the Fermi velocity can be assumed to be much larger than the sound 
velocity. Then, the propagator Q^(k, r — r') is a fast changing function of time with respect 
to the boson propagators, so that in Eq. (JT5j) the static limit of the effective interaction 
can be taken. In this case, the boson-induced interaction is the space-Fourier transform 
( times —A 2 ) of the static correlation function for density fluctuations. In the following 
calculations, this second approximation will be adopted, since it is closer to the physical 
situations realized in experiments. Thus, in the frequency representation, the equation 
for the pairing field with P = 0, is given by the following expression 

A(k ' Un) ~ ~ J ^ + 4 2 1 + |A(k 1 ,^ ) p • (18) 

When evaluating the energy gap in the spectrum of elementary excitations, an analytical 
continuation of this equation to real frequencies should be made, in order to determine 
the real poles of the fermion Green function. However, for assessing the occurrence of a 
superconducting phase and for evaluating the order of magnitude of the energy gap, it is 
sufficient to take into account only the u — component of the pairing field, A(k). This 
quantity satisfies the equation 

dk 1 V e// (k-ki)A(ki) 



*< k > = -7(# -% + iA(M' • (19) 

Here V e ff(k — k x ) denotes the effective interaction at equal times, its explicit expression, 
with q = |k — ki|, is 

V eff (q) = -X 2 n° B [l + 2n' B (q) - ^[1 + 2n{q)\) - A 2 / ^ [n' B {k)[l + n' B (|k - q|)] + 



+ A l „ [1 + 2*(*)][1 + 2n(|k - q|)] , (20) 

where Ej. are the excitation energies of the Bose gas calculated within the Bogoliubov ap- 
proximation, E k = J [k 2 j (2m B )) 2 + 2^/n B k 2 / (2m B ), n' B (k) the occupation numbers of non- 
condensate bosons, and 

m = T^i ■ 

The first term on right hand side of Eq. (|2U|) is the contribution from states containing 
one virtual phonon, and it comes from the correlations of the products of the conden- 
date and noncondensate amplitudes taken at two different space points. The second term 
is the contribution from states containing two virtual phonons, and it is related to the 
noncondensate-density correlations. 

A few remarks should be made. First, according to the Pauli principle, the pairing field 
should be an odd function of k in order to change sign under the permutation of the two 
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correlated fermions. This excludes contributions to the pairing field from interactions acting 
only in the s-wave channel. So that, as previously remarked, the bare fermion-fermion 
interaction can be completely neglected in the present scheme. In addition, the function 
V e ff(q) has a simple pole, when q — > 0. However, this is a fictitious pole, since, if the 
the integrals in Eq. (|2(Jj) are replaced by a sum over the momenta, the boson states with 
vanishing momentum should be excluded. Actually, this restriction may be omitted, because 
it affects the integrand of Eq. (J 19)) only at a single point and the integration gives a finite 
result. Finally, it is worth noting that V e ff(q) does not represent a usual two-body potential, 
since an integration over the time is absent. 

In order to solve Eq. (fE?j) it is convenient to expand A(k) in spherical harmonics 

A ( k ) = J2 ]j^^^ A i,m(k)Y l m (n k ) , 

with only odd values of I contributing. To simplify calculations, an angle-average approxi- 
mation for |A(ki)| 2 can be made 

lACkOI 2 - A 2 ^) = T,^— 1 \ A Uh)\ 2 • (21) 

l,m 

With this approximation, the different m-components turn out to be uncoupled and all 
equal; in addition, they can be assumed to be real and positive, without loss of generality. 
The quantity A 2 (k) represents the averaged energy gap in the quasiparticle spectrum of the 
superconducting phase. 

The equations for the partial waves of the pairing field, 

can still be coupled via the term A 2 (ki). Here, the components Vi(k, k\) are defined as 

Vi(k,k 1 ) = J dxV eff (k 2 + k 2 1 -2kk 1 x)P l (x). 

They are bounded quantities for k = k\, in spite of the occurence of a pole when the 
argument of K//(<z) vanishes. 



III. RESULTS 

The temperature, density and relative concentrations of the gas mixture are assumed such 
that a relevant fraction of the Bose gas is in the condensate phase and the Fermi gas is still 
quite degenerate. To be more specific, we refer to the experimental environment described 
in Ref. [5|, where a mixture of 87 Rb (bosons) and 40 K (fermions) atoms, with a relevant 
ratio between the fermion and boson densities, rip/riB ~ 0.5, has been studied. For the 
homogeneous boson gas, a value T c = HOnK for the critical temperature of condensation 
has been chosen, in order to reproduce values of temperature and condensate fraction close 
to those reported in Ref. 0]. The corresponding boson density is ub = 1.45 x 10 13 cm~ 3 . 
This is a reasonable value for the average boson density obtained in the experiment of 
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FIG. 1: Induced effective interaction between fermions as a function of the momentum transfer 
q, for a mixture of 40 K and 87 Rb (solid line), r$p is the mean spacing between fermions. The 
contributions from states with one virtual phonon ( dashed line ) and with two virtual phonons 
( dot-dashed line ) are also shown. The condensate fraction is n B /n B = 0.6. 

Ref. ff|. Moreover, for the boson-boson and boson-fermion scattering lengths, the values 
&bb — 113ao and a B p = — 330ao are used, where ao is the Bohr radius. The range of the 
condensate fractions considered in the present calculations is 0.4 < n° B /n B < 0.8, and the 
corresponding interval of temperatures is 46nK < T < 87nK. With these values of the 
relevant parameters, the ratio between the Fermi velocity, defined as vp = y2 fip/mp, and 

the sound velocity in the Bose gas, c s = ^J'jn B /m B , turns out to be about 20. This value of 
v f/c s is sufficiently large to justify the approximation of Eq. (|20|) to the effective interaction. 
Moreover, the Fermi temperature corresponding to the density rip ~ 0.5ns, Tp ~ 320nK, is 
high enough with respect to the temperatures in the interval considered, so that the Pauli 
principle is still operating. 

The results for the effective interaction are discussed first. A simple inspection of Eq. 
fl20j) shows that the contributions from states with one virtual phonon and from states with 
two virtual phonons, approach the values —X 2 n B and —X 2 n' B , respectively, when q — > oo. 

In Fig. the effective interaction of Eq. ()2())1 is shown togheter with the separate con- 
tributions from one-phonon and two-phonon states, for n° B /n B = 0.6. We can see that the 
two contributions are of the same order of magnitude for q > 0.5 r^, with tqp denoting the 
mean spacing of the Fermi atoms, and that the asymptotic values are practically reached al- 
ready at q ~ 1.5 Tq F . This suggests that the contribution from states containing two virtual 
phonons cannot be neglected, unless the condensate depletion is very small. 

In Fig. |21 the effective interaction is displayed for three different values of the condensate 
fraction. From Fig. |2]we can see that, in the considered range of n° B /n B , for a fixed value of 
the boson density the magnitude of the effective interaction increases with the condensate 
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FIG. 2: Induced effective interaction between fermions as a function of the momentum transfer q, 
for three different values of the condensate fraction. The Bose-Fermi mixture is the same as in 
Fig. n Dashed line: n B /nB = 0.8; solid line: n° B /riB = 0.6; dot-dashed line: n° B /riB = 0.4. 

depletion, or, equivalently, with the temperature. This is due to the fact that the increase 
with the noncondensate fraction of the two-phonon contribution, overcomes the decrease 
with n° B of the one-phonon contribution. However, for higher values of the condensate 
depletion, the increase of V e ff(q) results considerably quenched. 

Concerning the solutions of the gap equation (}2*2*j) . it should first be remarked that the fast 
approach of the effective interaction to an almost constant value makes the odd components 
Vi(k,ki) vanish rapidly at high momenta. As a consequence, for evaluating the integral 
of Eq. (|22|) it is not necessary to introduce any cut-off on high values of k\. The explicit 
solutions of Eq. (}2~2~|) display two general features for the physical situations considered: they 
give, for all values of I, a small energy gap, with respect to the effective chemical potential 
of the of the fermions, Ai(k)/fi* F ~ 10~ 2 , and they give a non vanishing energy gap only 
for a single value of I at a time, without any coupling to the remaining components. The 
latter feature can be better understood by observing that, for small values of Ai(k), the 
most relevant contribution to the integral of Eq. ()22|) comes from a small domain about the 
effective Fermi momentum, k* F = ^2fi* F m F . This allows us to approximate all the quantities 
in the integrand, apart from by the values that they take at k\ = k* F . Thus, Eq. (122}) 
turns into the following closed equation for each value of k 
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FIG. 3: Energy gap for different /-channels, as a function of the wave vector k measured in units 
of kp = ^J2/j,* f nip. Solid line: I = 1; dashed line: I = 3; dot-dashed line: I = 5. The condensate 
fraction is n B /riB = 0.6, corresponding to a temperature of 70nK. The Bose-Fermi mixture is the 
same as in Fig. ^ 

It is easy to see that Eq. ()23j) can have nontrivial solutions only when Ai> = for I' ^ I. 
In addition, it is worth noting that the energy gaps given by Eq. ()23j) represent a very 
satisfactory approximation to the numerical solutions of Eq. (J22j) . 

The solutions of Eq. (|2*2|l indicate that transitions to different superconducting phases 
could occur. Each phase is characterized by a nonisotropic pairing field, whose angular 
dependence is given by a specific spherical harmonic. The corresponding energy gap contains 
only one term Ai(k). 

In Fig. El the energy gap is displayed for three different values of I at a fixed condensate 
fraction ( or temperature ). The fast decrease of the energy gap with increasing / corresponds 
to the analogous behavior of the effective interaction. While the peak shown by Ai(k) 
about k = k* F is due to the sharp increase of the effective interaction when the transferred 
momentum appoaches zero. This can be directly seen from the approximate equation (|23|) . 

The behavior with the condensate fraction, of the energy gap for I — 1 is reported in 
Fig. |3J All /-channels show a similar behavior. We see that, when n° B decreases, or, equiv- 
alently, the temperature increases, the gap Ai(fc) increases for all values of k. This feature 
essentially reflects the analogous behavior of the effective-interaction strength. However, the 
fermion gas becomes less and less degenerate when the temperature increases, so that, the 
corresponding decrease of the effective Fermi momentum tends to counterbalance the effect 
of the effective interaction. Actually, explicit calculations show that, for T ~ lOOnK corre- 
sponding to n° B /riB ~ 0.2, the gap Ai(fc) reverses its initial behavior. It should be observed 
that the approximation scheme proposed in this paper does not keep its validity for higher 
temperatures. In any case, the results of the present approach show that a superconducting 
phase in the p-wave channel can be achieved at temperatures > lOOnK. In Ref. a critical 
temperature of ~ 10~ 2 nK was indicated for p-wave pairing in a mixture of 6 Li and 87 Rb. 
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FIG. 4: Energy gap in the I = 1 channel for different values of the condensate fraction. Solid 
line: n° B /n B = 0.4 ( T = 87nK ); dashed line: n° B /n B = 0.6 ( T = 70nK ); dot-dashed line: 
n %/ n B = 0.8 ( T = 46nK ). The Bose-Fermi mixture is the same as in Fig. ^ 

For this mixture calculations performed within the present scheme and with the same values 
of parameters as in Ref. show that the behavior of the energy gap qualitatively does 
not differ from that of Fig. HJ The deep discrepancy between the results of the present work 
and those of Ref. [H should be ascribed to the fact that in the two approaches two opposite 
approximations have been used for the time dependence of the effective interaction. These 
approximations have already been discussed in Sec. IIB. 

IV. SUMMARY AND CONCLUSIONS 

In a mixture of ultracold Fermi and Bose gases, an interaction between fermions is induced 
by the exchange of elementary excitations of the Bose gas. Due to this supplementary 
interaction, a superconducting phase can appear in the Fermi component of the mixture. 
The mechanism is analogous to that accounting for superconductivity in ordinary metals. 
In this paper, a procedure to determine the induced interaction and to assess the possibility 
of a phase transition has been proposed. For simplicity, a spatially homogeneous system 
has been considered. Both the gases have essentially been treated within a mean field 
approximation appropriate for the physical conditions realized in recent experiments: the 
values of temperature and densities are such that the Bogoliubov approximation for the 
Bose gas can be justified and the Fermi gas can be considered degenerate. Moreover, since 
the experiments concern gases confined by means of magnetic traps, in the present approach 
the atoms of both species have been considered completely polarized. This circumstance 
has allowed the bare fermion-fermion interaction to be neglected. 

The same formalism, based on functional methods, has been used in all stages of the cal- 
culation. This allows us to introduce the various approximations to the effective interaction 
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between fermions in a self-consistent way. In particular, the inclusion of the contributions 
from states containing two virtual phonons and the static approximation to its time depen- 
dence. This is more appropriate for the physical situations of actual experiments than the 
usual instantaneous approximation. Both terms of the interaction, with one and two virtual 
phonons, are attractive and of the same order of magnitude. They can give rise to the pair- 
ing of correlated couples of fermions with different values of the center of mass momentum. 
In this paper, only the pairing field with vanishing center of mass momentum has been 
considered, thus preserving the invariance of the system under space translations. Since 
the paired fermions lie in the same magnetic state, the pairing field is constrained to be an 
odd function of their relative momentum. Then, it contains only odd-Z components, when 
expanded in spherical harmonics. A particular feature of the corresponding gap equations is 
that they show solutions with only one non vanishing partial wave at a time. Then, differ- 
ent nonisotropic phases may occur in the Fermi component of the mixture, with the pairing 
fields decreasing their strength when their multipolarity increases. It is worth noting that 
the energy gaps increase with the temperature until a value of T ~ lOOnK, corresponding 
to a high depletion of the Bose condensate n° B /n B ~ 0.2. The approximations introduced 
in the present approach become less reliable for higher values of the temperature, for the 
values of boson and fermion densities assumed here. Nevertheless, we can expect that a 
superconducting phase can occur in a Bose-Fermi mixture at T > 100 nK. This value is a 
relevant fraction of the Fermi temperature. 
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